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!h ■ ABSTRACT 
^— > . 

VI . 

. The morphological and spectroscopic properties of the giant (>60 kpc) Lya neb- 

ulae associated with three radio galaxies at z -2.5 ( MRC 1558-003, MRC 2025-218 
\ and MRC 0140-257) have been investigated using integral field spectroscopic data 

obtained with VIMOS on VLT. 
\j£ 1 The morphologies are varied. The nebula of one source has a centrally peaked, 

H . rounded appearance. In the other two objects, it consists of two spatial components. 

■ The three nebulae are aligned with the radio axis within <30°. The total Lya lumi- 
nosities are in the range (0.3-3.4) x 10 44 erg s -1 . The Lya spectral profile shows strong 

Tij-j- . variation through the nebulae, with FWHM values in the range —400-1500 km s _1 

■ and velocity shifts V ffset —120-600 km s _1 . 

£^ \ We present an infall model that can explain successfully the morphology, size, 

. surface brightness distribution and the velocity field of the Lya nebula associated with 

MRC 1558-003. It can also explain why Lya is redshifted relative to other emission 
. t-h ! lines and the FWHM values of the non resonant Hell line. This adds further support 

■ to our previous conclusion that the quiescent giant nebulae associated with this and 
^ \ other high redshift powerful radio galaxies are in infall. A problem for this model is 

the difficulty to reproduce the large Lya FWHM values, which might be consequence 
of a different mechanism. 

We have discovered a giant (—85 kpc) Lya nebula associated with the radio galaxy 
MRC 0140-257 at z =2.64. It shows strikingly relaxed kinematics (FWHM<300 km 
s _1 and Voffset ^$120 km s _1 ), unique among high z (>2) radio galaxies. 

Key words: galaxies: active; galaxies: high redshift; galaxies: individual: MRC 1558- 
003, MRC 2025-218, MRC 0140-257 



1 INTRODUCTION 

Powerful high z radio galaxies [z >2, HzRG) are often sur- 
rounded by giant Lya nebulae whic h can extend for more 
than 100 kpc (e.g. McCarthy et al. [l990a; Reuland et al. 
2003, Villar-Martm et al. 2003) and som etimes beyo nd the 
radio structures (e.g. Eales et al. [j.993; Kurk et al. 120021 : 



* Based on observations carried out at the European South- 
ern Observatory, Paranal (Chile). Programs 075.B-0212(A) and 
073.B-0189(A). 



Maxfield et al. l2002h . These are clumpy, irregular (with 
features such as fil aments, plumes, ionization cones, e.g. 
Reuland et al. l2003r) a nd often aligned with the radio axis 
(McCarthy et al. 1995). They are characterized by extreme 
kinematics, with measure d FWH M >1000 km s _1 (e. ^Mc- 
Carthy, Baum & Spinrad [l996|; Villar-Martm et al. I2003T ). 
compared with values of a few hundred in lo w-re dshift ra- 
dio galaxies (e.g. Tadh unter, F osbury & Quinn|l989|; Baum, 
Heckman & van Breuge l 19901 ) . There is strong evidence that 
interactions between the radio structures and the ambient 
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gas, produce an outflow responsible for such extreme kine- 
matics (e.g. Humphrey et al. 120061 van Ojik et al. 1997). 

The nebulae have typical values of ionized gas mass 
^10 9-10 M , Lya luminosities ^severalx 10 43-44 erg s _1 
and densit ies n e ~ few to several hund red cm -3 (e.g. 
McCarthy [l993l . Villar-Martin et al. l2003h . They emit a 
rich emission line spectrum dominated in the optical (UV 
rest frame) by Lya followed by CIVA1550, HeIIA1640 and 
CIII]A1909 (CIV, Hell and CIII] hereafter). Such a line spec- 
trum reveals high levels of metal enrichment and excitation 
mechanisms mostly related to the nuclear activity, at least 
in th e direction along the radio structures (e.g. Vernet et al. 
l200lh . 

In addition to the highly perturbed gas, HzRG are often 
embedded in giant (often ^100 kpc), low surface brightness 
nebulae of metal rich, ionized gas with q uiesce nt kinematics 
(Villar-Martm et al. l200l Reuland et al. 120071 ). i.e., not per- 
turbed by interactions with the radio structures. We have 
recently shown that these quiescent nebula e are inf ailing to- 
wards the central region (Humphrey et al. l2007h . 

There are only several narrow band Lya images of 
HzRG and most spectroscopic studies have been performed 
with the long-slit technique, with the slit aligned with the 
radio structures. These studies have been seriously limited 
by, respectively, the lack of spectral information, or the lack 
of spatial information in directions other than the radio 
axis. For this reason, we are carrying out an observational 
program of 3D integral field spectroscopy of powerful radio 
galaxies at z -2-3 with VIMOS on VLT and PMAS/PPAK 
on the 3.5m telescope in Calar Alto Observatory. The main 
goal is to characterize the morphological, kinematic and ion- 
ization properties of the extended ionized gas in two spatial 
dimensions. 

In this paper, we present results obtained for MRC 
1558-003 (*=2.53), MRC 2025-218 (z=2.63) and MRC 0140- 
257 (£=2.64), based on VIMOS- VLT d ata. Results on MRC 
2104-242 (2;=2.49) can be found in IVillar-Martm et al. I 
<|2006h . A similar study of 4C40.36 (z=2.27) and 4C48.48 
(2=2.34) based on PMAS/PPAK data (3.5m telescope, 
Calar Alto Observatory) will be presented in Sanchez et al. 
(2007, in prep.). 

A Q A = 0.73, Q m 0.27 and H 62 km s" 1 Mpc" 1 
cosmology is adopted in this paper (Sandage et al. l2006h . 



2 OBSERVATIONS AND DATA REDUCTION 

The observations (program 075.B-0212(A)) were made on 
UT 2005 July 28, 29 and 30 using the Visible MultiObject 
Spectrograph (VIMOS, Le Fevre et al. l2003h , on the Nas- 
myth focus of the UT3 VLT. The instrument is equipped 
with an integral field unit with 6400 microlenses coupled to 
fibres. For the configuration selected by us, the number of 
fibers in use is 1600, covering 27" x 27" on the sky with a 
0.67" sampling. The HRbiue grating was used, with an effec- 
tive wavelength range ^4150-6200 A, an instrumental profile 
of FWHM 1.7±0.2 A and a pixel scale of 0.5 A pixel" 1 . 

The exposure time on each target was 7.3 h (22 x 1200 s) 
on MRC 1558-003, 10 h (30x1200 s) on MRC 2025-218 and 
8.3 h (25x1200 s) on MRC 0140-257. In addition, we had 
4.5 h (9x1800 s) on MRC 1558-003 obtained in J une 2 004 
(program 073.B-0189(A), see Villar-Martm et al. 120061 for 



a description of these observations). The total integration 
time on this source was therefore 11.8 h. 

The seeing full width at half maximum during the ob- 
servations was in the range ^0.44-1.32" (1st night), 0.4-1.2" 
2nd night) and 0.5" -3.0" (3rd night) for the 2005 observa- 
tions, and 1.0-1.4" for the 2004 data on MRC 1558-003. 

For each galaxy a dithering pattern was applied, with a 
maximum offset of ~3", and a range of dithering pointings 
between 3 and 7. 

The data were reduced using R3D (Sanchez 2006) and 
IRAF routines. The data were bias subtracted. The loca- 
tions of the spectra were traced on a continuum- lamp expo- 
sure obtained before each target exposure. The correspond- 
ing spectrum was then extracted for each fiber by coadding 
the flux intensity within a 5 pixel aperture, along the spec- 
tral dispersion axis, centered on the estimated location of 
the fiber centroid. 

The wavelength calibration was performed using arc 
lamp spectra and the telluric emission lines in the science 
data. The wavelength solution was stable within a fraction 
of 1 spectral pixel (0.5 A) across the whole spectral range 
and the whole field of view. 

The wavelength calibration in the blue (i.e. the Lya 
spectral region) was problematic due to the lack of bright 
emission lines in the arc lamp and sky spectra. Residual, 
artificial shifts in A of up to 6 A were measured for Lya rel- 
ative to CIV and Hell in MRC 1558-003 and MRC 2025-218. 
We have been able to correct for this effect by comparing 
with available Keck long- s lit spectroscopy (Villar-Martm et 
al. 120031 Humphrey I 20041 ). to an accuracy of <1 A. In any 
case, such uncertainty will not affect the results presented 
here. 

The fibre-to-fibre response at each wavelength was de- 
termined from a continuum-lamp exposure. 

After these basic reduction steps, a data cube was cre- 
ated for each exposure. The cubes for MRC 2025-218 were 
then recentered spatially at each wavelength by determin- 
ing the centroid of a nearby star in the VIMOS field of view. 
This recentering corrects for differential atmospheric refrac- 
tion. 

For the other two objects, there are no stars in the VI- 
MOS field of view. We used the information on the spatial 
shifts contained in the image headers to apply the corre- 
sponding shifts. However, we found that this technique left 
substantial spatial offsets. For this reason, we used the Lya 
peak of emission to recenter all cubes. Although this recen- 
tering is valid in the blue, it is not clear that it is also valid 
in the red. In particular, in the case of MRC 1558-003, we 
found a residual spatial shift of ~1 spaxel in the direction of 
the radio structures between the blue (Lya) and red (CIV, 
Hell) parts of the spectrum, which is much larger than found 
in previo us works (Villar-Martm et al. 1999; Humphrey et 
al. l2007h . For this reason, an additional correction was ap- 
plied. Uncertainties remain regarding the spatial centering 
in the direction perpendicular to the radio axis, although 
this is likely to be <1 spaxel or 0.67". This will not affect 
seriously the results presented here. 

For MRC 0140-257, no useful previous works or addi- 
tional data were available. Although the accuracy of the 
spatial centering between the red and the blue ends of the 
spectrum is uncertain, this will not affect the results pre- 
sented here, due to the nature of our analysis and the clear 
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separation between the spatial components of the Lya neb- 
ula. 

The cubes were then combined using IRAF tasks, mask- 
ing the broken and/or low sensitivity fibres. A 3a clipping 
algorithm removed cosmic rays. The sky background was 
estimated before subtraction by selecting the spectra of ob- 
ject free areas and c reating an interpolated datacube (using 
E3D, Sanchez 1 2 004 ). A spatial median smoothing using a 
4x4 spaxel box was applied. 

The nights were non-photometric. Observations of stan- 
dard stars were used to perform a relative calibration from 
blue to red. The absolute flux calibration for MRC 1558-003 
and MRC 2025-218 was done using available Keck long-slit 
spectroscopy along the radio axis. We extracted a pseudo- 
slit from the Vimos data as similar as possible to the slit 
aperture of the Keck spectra. The flux in different spectral 
windows was then compared and scaled. The final flux cali- 
bration has an accuracy of -20%. For MRC 0140-257, this 
could not be done, due to the unavailability of a high quality 
long-slit spectrum, but as we explain in §5.3, the agreement 
between the measured Lya flux and published values is rea- 
sonably good (—20%). 

Cross-talk effects are estimated to be negligible. The 
fibre to fibre contamination is expected to be <5% for adja- 
cent spectra in the CCD, dropping to less than 1% for the 
3rd adjacent spectra. 

To ove rlay the radio maps (Carilli et al. 1997, Pentericci 
et al. l2000h on the Lya images, a different method was used 
for each object depending on the information available in 
the data. For MRC 2025-218, we placed the radio core at 
the position of the continuum centroid in the VIMOS data. 
Since this is a broad line object (§3.2), it is reasonable to 
expect that this marks the location of the active nucleus 
(AGN). This is further supported by the fact that the UV 
continuum has an spatially unresolved component (§3.2). 

For MRC 1558-003, we positioned the radio core at the 
spaxel with the maximum Lya flux, which is expected to be 
shifted by a few tenths of a spaxel relative to the contin- 
uum centroid (Humphrey et al. 120071 ). This method would 
be incorrect if the AGN is spatially shifted relative to the 
continuum and/or line centroids. However, we do not expect 
this to be the case since this is a broad line object and more- 
over, it shows little evidence for Lya absorption (see below). 
The la uncertainty in the Lya-radio registration is —0.3" 
in both cases. 

For MRC 0140-257, we assumed that the radio core is 
placed between the two Lya spatial components, although 
it is not necessarily the case. This is the main source of 
uncertainty since the radio core could be shifted by 2" if it is 
associated with one of the two Lya blobs (§4.3). The impact 
of this assumption on the interpretation of our results will 
be discussed when relevant. 

2.1 Vimos sensitivity 

The main scientific goal of our VIMOS- VLT observational 
program is to study the properties of the giant nebulae as- 
sociated with a sample of HzRG. By isolating spectroscopi- 
cally and/or spatially the emission from the perturbed and 
the q uiescent gaseous components (e.g. Villar-Martm et al. 
2003), one of our priorities is to characterize the morpho- 
logical, kinematic and ionization properties of the quiescent 



gas. In this way we can study the gas without the strong 
distortions that the radio structures can imprint on its ob- 
served properties. Such study has the potential to provide 
critical information on the star formation and chemical en- 
richment histories a s well as the galaxy fo rmation pro cess 
(Villar-Martm et al. l200l Humphrey et al. I2006L l2007h . 

In order to achieve these goals, it is critical to detect 
the main UV lines (very especially Lya and Hell) with high 
signal/noise in the faint, outer regions of the objects, which 
are usually characterized by quiescent kinematics. Moreover, 
this would allow to study in detail the high surface bright- 
ness regions (at least several times brighter). Although usu- 
ally distorted by the radio structures, it should be possible 
to isolate spectroscopically the emission from the quiescent 
and the perturbed gas using especially the non resonant Hell 
line, as Villar-Martm et al. (|20Q3h . 

The faintest regions of our interest have often Lya sur- 
face brightness levels which are, within the errors, consis- 
tent with or even below the detection limit or our Vimos 
data (3a — 10 -17 erg s _1 cm -2 arcsec -2 ) for 8 -10 h ours 
exposure time. For comparison, van Ojik et al. (1996) de- 
tected in 4 hours the giant, quiescent, low surface bright- 
ness nebula (— 10 -17 erg s _1 cm -2 arcsec -2 ) associated with 
a radio galaxy at z =3.6 using long slit spectroscopic data 
obtained with EMMI (ESO Multi-Mode Instrument) on the 
3.5m NTT telescope (La Silla). With a 2.5" wide slit, the 
spectral resolution was of 2.8 A, not very different to that 
of our data. In the VIMOS data, not even Lya could be 
definitively detected from the faintest regions of our inter- 
est, which at least in MRC 1558-003 and MRC 2025-218 
we knew exist (§3). Hell and CIV are detected only from 
the highest surface brightness regions, but with not enough 
signal/noise ratio to perform an adequate kinematic and ion- 
ization analysis in two spatial dimensions. 

Thus, the investigation of the properties of the quies- 
cent gas has been possible only for those objects where there 
is no signature of interactions between the radio structures 
and the gas (MRC 1558-003 and MRC 0140-257, see §4) or 
objects with no spectroscopic blend between the perturbed 
an d the quiescent nebula (MRC 2104-242, Villar-Martm et 
al. l2006r ) . Moreover, in all three cases, Lya is not heavily ab- 
sorbed. For objecs where the emissions from the perturbed 
and the quiescent gas are blended and/or Lya is absorbed 
(i.e. an important fraction of HzRG), a more sensitive in- 
strument/telescope combination is needed. 



3 PREVIOUS RESULTS 

The three radio galaxies discussed in this paper, MRC 1558- 
003, MRC 2025-218, MRC 0140-257, belong to the Mo longlo 
Catalogue of high redshift radio galaxies (Large et al. 1981, 
McCarthy et al. [l990b). The objects were selected based 
on previous evidence for strong Lya emission. The two first 
objects were previously known to be associated with giant 
(>60 kpc) Lya nebulae (see below). 

3.1 MRC 1558-003 (z=2.53) 

This object is associated with radio struct ures which extend 
for -9" or 84 kpc (Pentericci et al. 2000) along a position 
angle (PA) of 75°. 
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Previous long-slit spectroscopic studies have shown that 
Lya is extended for at least ~14" (^130 kpc) along the slit 
with PA =7 2° we ll beyond the radio structures (e.g. Villar- 
Martm et al. 2003). CIV and Hell are extended across ^6" , 
while NV, which is fainter, has a more compact appearance. 
The optical continuum is also extended. 

ISAAC near infrared spectroscopy revealed very 
broad Ha (FWHM~11700±900 km s~\ Humphre£2Q0l; 
Humphrey et al. 2007, in prep.), evidence for an obscured 
broad line region (BLR). According to the unified scheme of 
quasars and radio galaxies (Barthel 1989), the detection of 
BLR emission implies that the ionization cones axis is at a 
small angle to line of sight (<45°). This is further suggested 
by the detection of a clear one-sided radio jet (see Fig. 1). 

High resolution optical spectroscopy (FWHM^1.7 A) 
showed no absorption features in the Lya spect ral profile 
integrated along PA=85° (van Ojik et al. Il997h . Absorp- 
tion features w ere no t found either along PA=75° by Villar- 
Martin et al. (|2003T ). although the spectral resolution was 
rather low in this case (FWHM-11 A). The small impact 
of absorption effects is also suggested by the la rge Ly a /Hell 
values measu red ac ross the nebula (Humphrey [ 2004 Villar- 
Martm et al. I2007T ) 

For this object no HST or emission line images are avail- 
able in the literature. 



3.2 MRC 2025-218 (z=2.63) 

This radio galaxy is associa ted w ith a small radio source 
(-5" or 46 kpc, Carilli et al. Il997h . 

There is strong evidence for an obscured BLR in MRC 
2025-218. The morphology of the object at optical and near 
infrared wavelengths is dominated by a point source (Penter- 
icci et al. 1999, 2001; McCarthy, Person & West 1992) . Broad 
CIV and broad Ha (FWHM-6000 kins" 1 ) revealed by op- 
tic al (Vi llar-Martin et al.ll999h and near infrared (Larkin et 
al. l2000l : Humphre^ T 20041 : Humphrey et al. 2007, in prep.) 
spectra confirm this interpretation. The UV rest frame con- 
tinuum emission shows also an extended (^6") diffuse com- 
ponent, which is well aligned with the radio axis (Pentericci 
et al. 1999). 

The Lya emission is distributed bimodally (McCarthy 
et al. Il990bl ) north-south with two peaks that correspond 
roughly to the two radio lobes. The authors measured a to- 
tal extension of ^4" or 37 kpc. Keck long-slit optical spec- 
troscopy (Humphre}|_200J) reveals that Lya is extended up 
to ~9" or 83 kpc. CIV, Hell and CIII] are also spatially 
extended. CIV and CIII] each have a spatially unresolved 
component centered at the continuum centroid and also an 
underlying extended component spanning ^5" along the ra- 
dio axis. NV is spatially u nresol ved. 

Villar -Martin et al. ] ([1999) detected absorption in 
the spectrum of MRC 2025-218 for CIV, CIIA1335, 
SiIVAA1393.8,1402.8 and, maybe, OIA1302.2+ SiIIA1402.8. 
The steep Lya profile on the blue side of the line was also 
proposed to be a signature of absorption. 

3.3 MRC 0140-257 (z=2.64) 

This object is associated w ith a small double radio source 
(4.2" or 39 kpc, Carilli et al. ll997l ). In the optical, it appears 



as a faint galaxy spatially extended and aligned with the 
radio source (McCarthy, Persson & West 1992). 

Lya and Ha are the two lines d etecte d in previous spec- 
troscopic studies (McCarthy et al. 1991; Eales & Rawlings 
1996). It has not been reported whether the lines are spa- 
tially extended. 

The F160W filter NICMOS-HST image (rest frame 
spectral window 3850-4950 A, Pentericci et al. l200ll ) shows 
two peaks of emission with almost the same flux, closely 
aligned with the radio axis. The authors suggest that this 
morphology could be due to a dust lane. There are a few 
fainte r clumps within 2 -3" . 

lEales Rawlingj (|l996h report Lya/Ha=0.31, which 
is ^30 times lower than standard case B photoionization 
model predictions. As the authors explain, this value is very 
uncertain, since they had to assume a fixed value of the 
[Nil] /Ra ratio (0.5) and no aperture corrections to the line 
fluxes were attempted. In spite of these uncertainties, the 
very large discrepancy with the model predictions suggest 
that Lya absorption/extinction is present. 

For this object no emission line images are available in 
the literature. 



4 RESULTS 

We present below the most relevant results obtained for the 
three radio galaxies investigated here. 

4.1 MRC 1558-003 (z=2.53) 

The Lya nebula 

The VIMOS Lya+continuum image of MRC 1558-003 was 
created by adding the monochromatic images of the object 
within the [4285-4315] A range. A continuum image was sub- 
tracted to keep the line emission only. This was created by 
combining the images extracted from two adjacent spectral 
windows at both sides of Lya, with the same spectral width 
as the line+continuum image. The resulting Lya image is 
shown in Fig. 1 (see also Fig. 2), with the 8.2GHz radio 
contours overlaid (Pentericci et al. l2000h . 

The nebula extends across ~9"x7.5" or 84x70 kpc 2 
(but see §3.1). It is characterized by a strongly peaked 
rounded morphology. Emission from the nebula is detected 
only within the radio structures. The faintest Lya emis- 
sion we detect on the VIMOS image has surface brightness 
~10 -17 erg s _1 cm -2 arcsec -2 (3a values). 

The angle between the nebular axis and the radio struc- 
tures was measured using the position angle on the sky of 
the longest dimension of the Lya nebula and the position 
angle of the line between the two brightest radio hot spots. 
The nebula is misaligned by ^30±5° relative to the radio 
axis. 

As we found for MRC 2104-242 (Villar-Martin et al. 
2006), Lya emission is detected outside any pl ausibl e ioniza- 
tion cones with opening angle ^90° (Barthel 1989). Seeing 
effects are not likely to be responsible for this Lya emission 
since the observations of this object were carried out under 
seeing conditions of FWHM< 1.5". Since this is a broad line 
object, it is possible that the broad rounded morphology of 
the Lya nebula is due to orientation effects (see §5.3). If the 
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CIV + Lya contours 



-8 -6 -4 -2 2 4 6 8 

Offset in RA (") 

Figure 1. Lya nebula (grey scale) associated with MRC 1558-003 
with the 8.2GHz radio contours overlaid. The maximum extension 
of the nebula, as measured from the VIMOS data, is ~9"x7.5" 
(84x70 kpc 2 ). The nebula is misaligned by ~30±5° relatively to 
the radio structures. North-up; East-left. 



ionization cone axis lies close to the line of sight, a broader, 
rounder morphology is expected. 

The total Lya flux integrated over the nebula is 



4.8xlO~ 10 erg s" 
3.4xl0 44 erg s _1 . 



cm corresponding to a luminosity of 



CIV and Hell morphologies 

The CIV and Hell images were created by adding the 
monochromatic images of the object within the spectral win- 
dows [5455-5485] A and [5780-5800] A respectively. Contin- 
uum images adjacent in wavelength to each line and of the 
same spectral width were subtracted to keep the line emis- 
sion only. 

The resulting images are shown in Fig. 2. The Lya con- 
tours have been overplotted on the top panels for compar- 
ison. The CIV line is spatially extended with a maximum 
projected size of —7.5" . The outer, low surface brightness 
regions extend in the same direction as Lya. Hell, which is 
fainter, appears more rounded and compact, although it is 
spatially resolved with a maximum extension of —3". 

The inner (°°2 // ) regions of the CIV nebula seem to 
be extended along an axis (roughly E-W) which is shifted 
by — 15° ant i- clockwise relative to the axis denned by the 
outer regions and aligned within a few degrees with the ra- 
dio structures. The Hell emission seems to be extended also 
in this direction. There is some hint that this could also be 
the case for the inner Lya nebula. However, higher spatial 
resolution data would be necessary to investigate whether 
there is a real rotation of the nebular axis as we move out- 
wards, which could be a consequence of an inhomogeneous 
distribution of material. 

Continuum is also detected, although the image is too 
noisy to characterize its morphology and determine its spa- 
tial centroid. 

Spatially extended emission line spectroscopy 

We show in Fig. 3 the spatial maps of the FWHM (cor- 
rected for instrumental profile) and velocity shift Voffset of 
the Lya line. These values have been measured fitting Gaus- 
sian profiles spaxel to spaxel. V f fset has been computed rel- 
ative to the Lya emission at the spatial line centroid. The 




E <~ 



CIV 




13.4" (~ 125 kpc) 



Hell + Lya contours 




Hell 



Figure 2. MRC 1558-003: Comparison between the CIV, Hell 
and Lya: spatial distributions (see electronic manuscript for the 
colour version fo the figures). The underlying continuum has been 
subtracted from all images. Top panels: The Lyo; morphology is 
shown with contours overlaid on the CIV (top-left) and Hell (top- 
right) images (colour scale). The bottom panels show the CIV and 
Hell morphologies. Lya contours: (0.15, 0.6, 1.0, 1.6, 2.0, 3.2, 4.4, 
4.4, 5.0) xlO -16 erg s _1 cm -2 arcsec" 
4.1, 5.5, 7.0, 8.3)xl0" 
(1.8, 2.4, 3.2, 3.9, 4.6, 5.3)xlQ- 17 
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errors on the FWHM and V ff se t are estimated to be, in 
general, <100 km s _1 and <40 km s _1 respectively. 

Although the line profile is asymmetric at some spatial 
positions, this method allows a first order study of the spa- 
tial variation of the Lya spectral profile across the nebula. 
The line spectral profile shows strong spatial variations, as 
is obvious from Fig. 3. 

The first striking characteristic is that Lya presents 
an increasing redshift towards the nebular centroid, with 
a maximum shift in velocity of —400 km s _1 at this position 
relative to the outer regions. 

The FWHM varies between -450 and 1350 km s" 1 
across the nebula. Three distinct regions (A, B and C in 
Fig. 3) can be isolated in the FWHM map according to the 
line width. The maximum values are measured in region 
B (in the range 1050-1340 km s _1 ), which runs approxi- 
mately along the radio axis and contains the nebular cen- 
troid. The Hell FWHM measured from the integrated spec- 
trum is 600±100 km s _1 . Region A shows narrower Lya, 
although the line is still quite broad (FWHM in the range 
700-1030 km s" 1 ) while the Hell FWHM is 650=b200 km 
s _1 in the integrated spectrum. At some positions Hell is as 
narrow as 450±50 km s _1 (consistent with Villar-Martm et 
al. I2003T ) . Lya trends to be notably br oader than Hell (this 
was also found by Villar-Martm et al. 120031 along the radio 
axis) . 

Region C is clearly different both in line width 
(FWHM(Lya)=650±50 km s~\ Fig. 4, bottom) and veloc- 
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Figure 3. MRC 1558-003: 2-dim Lya spectral properties (see 
electronic version for colour figures). The FWHM (left) (corrected 
for instrumental broadening) and the velocity shift (right) relative 
to the line emission at the Lya spatial centroid are shown. Values 
in km s _1 . Lya: flux contours are overplotted. Regions A, B and 
C (see text) are identified. (Notice that the spaxels on which the 
FWHM and V Q / f se t are determined are not the same and depend 
on the signal to noise ratio). 

ity shift. This region shows the largest blueshift relative to 
the Lya centroid (350±20 km s _1 for the integrated spec- 
trum). CIV is detected in region C, although noisy, and it 
has FWHM=600±100 km s _1 . 

Lya absorption 

We do not find definitive Lya absorption features in 
the VIMOS data across the nebula. This is consistent with 
previous studies (§3.1). 



4.2 MRC 2025-218(z=2.63) 

The Lya nebula 

The VIMOS Lya+continuum image of MRC 2025-218 was 
created by adding the monochromatic images of the object 
within the [4385-4445] A range. A continuum image ex- 
tracted from an adjacent line- free region of the same spectral 
width was subtracted to keep the line emission only. The re- 
sulting Lya image is shown in Fig. 4, with the VLA 8.2GHz 
radio contours overlaid (Carilli et al. Il997l ; see also Fig. 5). 

The Lya maximum extension is ^6.0" or ^55 kpc (but 
see §3.2), very similar to the radio source size . The faintest 
Lya emission we detect on the VIMOS image has surface 
brightness ^1.5xl0 -17 erg s _1 cm -2 arcsec -2 (3a values). 
The maximum extension in the direction perpendicular to 
the radi o structures as measured from the VIMOS data is 
~4". As lMcCarthv et al. I ([1990b) already pointed out, the 
nebula is aligned with the radio structures within a few de- 
grees. The bimodal distribution described in §3.2 is clearly 
seen. 

The total Lya flux integrated over the nebula is 
2.1 xlO -15 erg s _1 cm -2 corresponding to a luminosity of 
1.7xl0 44 erg s _1 . 

CIV and continuum morphologies 

The CIV image (Hell is too faint) was created by adding 
the monochromatic images of the object within the spectral 
window [5610-5650] A and subtracting the adjacent contin- 
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Figure 4. Lya nebula (grey scale) associated with MRC 2025-218 
with 8.2GHz radio contours overlaid. The maximum extension of 
the nebula as measured from the VIMOS data is ~6" or 55 kpc. 
The nebula is similar in size to the radio source and is closely 
aligned with it (but see text). North-up; East-left. 



uum. The continuum image of the object was created by 
integrating across the [5670-5850] A spectral window. 

The resulting images are shown in Fig. 5. The Lya con- 
tours have been overplotted (top panels) for comparison. 

CIV extends for ^4" in the same direction as Lya (N- 
S). The bimodal distribution shown by Lya is distinguished 
in the CIV image. Interestingly, the continuum is spatially 
unresolved in this direction, but is barely resolved in the E- 
W direction, with a FWHM of ~1.1±0.1", compared with 
FWHM=1.5±0.2" for the star in the field. The continuum 
centroid is located between the two Lya spatial components, 
i.e., where this line presents a minimum. 

Spatially extended emission line spectroscopy 

We show in Fig. 6 the spatial maps of the FWHM and 
velocity shift V ff se t of the Lya line. As before, these values 
have been measured by fitting Gaussian profiles to the line 
spaxel by spaxel. V ff se t has been computed relative to the 
Lya emission at the position of the continuum centroid. The 
errors on the FWHM and V ff se t are estimated to be, in 
general, <100 km s _1 and <35 km s _1 respectively. 

This analysis reveals two clearly distinct regions (north 
and south), which are different both in FWHM and V ff se t- 
These are coincident with the two Lya spatial components 
discussed above. The continuum centroid (located between 
both regions) shows intermediate FWHM values. The Lya 
spectra extracted from the apertures N and S (north and 
south regions, Fig. 6) and the region in between (2 spaxels, 
see aperture C in Fig. 6) are shown in Fig. 7 (left pan- 
els). Both regions emit also CIV and Hell (Fig. 7, top-right 
panel). The CIV-Hell spectra (Fig. 7, top right panel) have 
been extracted using ony the highest flux spaxels in the N 
and S regions (5 and 4 spaxels respectively) in order to max- 
imize the signal to noise ratio for the lines. The differences 
in FWHM and V ff se t seen in Fig. 6 are obvious also here. 

The S region shows the narrowest Lya spectral pro- 
files (FWHM in the range 500-650 km s _1 ) and the gas is 
blueshifted relative to the continuum centroid (values in the 
range — 100 to -200 km s -1 ). CIV (which is a doublet and, 
therefore, expected to be intrinsically broader) and Hell are 
also rather narrow, with FWHM=650±30 and 500=b30 km 
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Figure 5. MRC 2025-218: Comparison between the CIV, contin- 
uum and Lya spatial distributions. Top panels: The Lya morphol- 
ogy is shown with contours overlaid on the CIV (continuum sub- 
tracted, left) and continuum (right) images. Bottom panels: CIV 
(left) and continuum (right) images. CIV is extended in the same 
direction as Lya. The continuum centroid is located between the 
two Lya spatial components, where the line flux presents a min- 
imum. Lya contours: (0.2, 0.5, 1.1, 1.4, 1.7, 2.1, 2.7, 3.0)xl0- 16 
erg s _1 cm -2 arcsec -2 . CIV contours: (1.4, 2.0, 2.5, 2.8, 3.2, 
3.5) xlO -17 erg s _1 cm -2 arcsec -2 . Continuum contours: (1.7, 
3.2, 4.5, 6.0, 7.5) xlO" 16 erg s 
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Figure 7. MRC 2025-218. Top: Comparison between the Lya 
(left) and CIV-Hell (right) spectra from the northern and south- 
ern regions (see text). The southern CIV-Hell spectrum has been 
shifted artificially in the flux scale for clarity. The two lines are de- 
tected in both regions. Bottom: Lya (left) and CIV (right) spectra 
at the position of the continuum centroid. Notice the very broad 
wings of both lines, signature of the broad line region, as well as 
the absorption features. In spite of the clear difference in FWHM, 
the shift in A of the Lya line and the spatial separation between 
the two regions, notice the exact coincidence in wavelength of 
the sharp edge in the blue wing of the line for the three spectra. 
Fluxes in units of 10 — 16 erg s — 1 cm -2 A -1 . 



N 




4.7" 



500 600 700 BOO 900 IOOO 1100 -300 -200 -lOO O lOO 200 300 



Figure 6. MRC 2025-218: 2-dim Lya spectral properties. The 
FWHM (left) (corrected for instrumental broadening) and the 
velocity shift (right) relative to the line emission at the continuum 
spatial centroid are shown. Values in km s _1 .The location of the 
continuum centroid is shown with an 'x'. Lya flux contours are 
overplotted. Regions N, S and C (continuum centroid) discussed 
in the text are overplotted. 

s _1 respectively, as measured from the spatially integrated 
spectrum. 

The N region is characterized by broader emission lines 
(FWHM in the range -850-1070 km s -1 for Lya) and the 
gas is redshifted relative to the continuum centroid (V ff se t 
in the range - +50 to +250 km s" 1 ). The CIV FWHM is 



1220+40 km s -1 . Hell is too noisy in the VIMOS spectrum 
to measure its FW HM. However, Keck long-slit spectroscopy 
fHumphrev I 20041 ) of this region implies FWHM= 1100+ 100 
km s _1 . Therefore, the three lines in the N region are very 
broad compared with the S region. 

Although absorption plays an important role in the Lya 
spectral differences across the nebula (see below), kinematics 
is also clearly having an effect, since the CIV and Hell lines 
show marked differences between the N and S regions. 

The Lya and CIV spectra at the position of the con- 
tinuum centroid are shown in Fig. 7 (bottom panels). Very 
broad wings (presumably from the BLR, see §3.2) are seen 
in the line profiles of both Lya and CIV. 

Lya and CIV absorption 

The Lya and CIV lines are clearly absorbed in this 
object (Fig. 7), as w as alr eady discussed by other authors 
(Villar-Martm et al. 1999). Since this paper is mostly fo- 
cused on the properties of the emission line nebulae, we will 
present here a general description and defer a more detailed 
analysis of the absorbers for another publication (Humphrey 
et al. 2007, in prep). 

Some peculiar properties of the Lya spectral profile 
are a consequence of absorption rather than kinematics: the 
multiple peaks, the fact that the flux drops below the contin- 
uum level at some wavelengths, the sharp, almost vertical 
edge of the blue wing of the line, and the identical wave- 
length over the nebula at which this sharp edge is measured 
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Figure 8. Lya nebula associated with MRC 0140-257 with 
8.2GHz radio contours overlaid. Two main Lya spatial compo- 
nents (Al and A2) are identified in the image very closely aligned 
with the radio structures. The size of the nebula is ~9 /x or 83 kpc. 
North-up; East-left. 



Continuum 




Continuum + Lya contours 


1 




■ 


— 


















i 








> 












> 



















13.4" (-124 kpc) 



Figure 9. MRC 0140-257: Comparison between the contin- 
uum and Lya spatial distributions. Left: continuum images. 
Right: With Lya contours overlaid. Lya contours: (0.2, 0.4, 
0.6, 1.0) xl0 — 16 erg s — 1 cm -2 arcsec -2 . Continuum contours: 

9 — 9 



(2.1,3.5,4.2,4.9) xl0 _lb ergs -1 cm -2 arcsec 



(see below) are most naturally explained by absorption. Sev- 
eral absorption features are detected in both Lya and CIV 
(Fig. 7). 

Lya absorption is detected in more than 30 spaxels. 
We set a lower limit to the size of the main absorber of 
-4.7 // x3.5 // or -43x32 kpc 2 . It is possible that this ab- 
sorber covers the Lya nebula completely. Since absorption 
is detected in CIV as well, this implies that the absorbing 
gas is highly ionized. 

A striking characteristic is that the sharp edge of the 
Lya blue wing (see Fig. 7) happens at almost exactly the 
same wavelength (4403. 3±0. 5 A) in all spaxels where we 
have been able to measure it (—20), revealing very little 
kinematic structure of the main absorber along the line of 
sight and accross its whole spatial extension. 

4.3 MRC 0140-257 (z=2.64) 

The Lya nebula 

The VIMOS Lya+continuum image of MRC 0140-257 was 
created by adding the monochromatic images of the object 
within the [4410-4430] A range. An adjacent continuum im- 
age was subtracted to keep the line emission only. The re- 
sulting Lya image is shown in Fig. 8 with the VLA 8.2GHz 
radio contours overlaid (Carilli et al. 1997). 

Two main spatial components (Al and A2 in the figure) 
are identified in the image, which are aligned with the ra- 
dio structures within a few degrees. We have assumed that 
the radio core is located between Al and A2 (but see §2). 
The Lya spatial centroids of Al and A2 are separated by 
— 3.5" or 32 kpc. The maximum extension of the nebula (Al 
and A2) is 9" or 83 kpc. It, therefore, extends well beyond 
the radio structures. This conclusion is independent of the 
location of the radio core along the A1-A2 line (§2). The 
faintest Lya emission we detect on the VIMOS data has 
surface brightness — 10 -17 erg s _1 cm -2 arcsec -2 (3a val- 
ues). 

The total Lya flux integrated over the nebula is 
3.6xl0 -16 erg s -1 cm -2 , which, taking aperture correc tions 
into account, is consistent with [McCarthy et al. I (|l99lh . It 
corresponds to a luminosity of 2.9 xlO 43 erg s -1 . 

Continuum morphology 



A continuum image was created by collapsing the data 
cube across the continuum spectral windows [5270-5550] A 
and [5640-5800] A, i.e. rest frame range -[1450-1530] A and 
[1555-1600] A respectively. For co mpari son, the optical im- 
age of McCarthy, Persson & West (|1992i ) was obtained using 
a Gunn-Thuan filter, covering the rest frame range [1675- 
1920] 0. We detect a faint source, which overlaps partially 
with component Al (Fig. 9). Due to the large separation in 
wavelength between the continuum and Lya images, we can- 
not discard that the spatial shift between the continuum 
centroid and Al is partially artificial (see §2). 

Although the source is apparently extended in the N- 
S direction, this is dubious. The faintest structures to the 
north are consistent within the errors with the image detec- 
tion limit. On the other hand, there are noise residuals at 
similar flux levels at different positions accross the VIMOS 
field of view. 

Spatially extended emission line spectroscopy 

We show in Fig. 10 the spatial maps of the FWHM 
and Voffset of the Lya line. V ff se t has been measured rel- 
ative to the Lya emission between Al and A2. The erros 
on Voffset are estimated to be <40 km s _1 . The errors on 
the FWHM are in the range 60-100 kins" 1 . The mam rea- 
son for these large relative errors is the uncertainty on the 
continuum level, due to the noise (this object is fainter) and 
the presence of an underlying broad component (see below) . 
The line is narrow, with FWHM<500 km s _1 accross the 
whole nebula, compared with typical values of HzRG. 

The Lya spectra integrated over the highest signal 
to noise spaxels of Al (16 spaxels) and A2 (14 spaxels, 
see Fig. 10, green lines) are shown in Fig. 11 (right pan- 
els). In both spatial components, the Lya spectral pro- 
file is dominated by a strikingly narrow component, with 
FWHM=250±50 and <120 km s" 1 for Al and A2 respec- 
tively. An underlying broad component seems to be also 
present in both components. This is most clearly detected 



1 The authors used a 2.5m telescope and an exposure time of 
2700 sec to obtain this image. 

2 Using a single Gaussian, as the fits used to produce Fig. 11 
(left) the narrow peak is broadened because of the broad wings 
and the derived FWHM have values of up to ~500 km s — 1 . 
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Figure 10. MRC 0140-257: 2-dim Lya spectral properties. The 
FWHM (left) (corrected for instrumental broadening) and the 
velocity shift (right) relative to the line emission at the inter- 
mediate spatial position between Al and A2 are shown. Values 
in km s _1 . Only coloured (i.e., not white) spaxels in the FWHM 
map have measured FWHM values. V j jset is shown for the same 
spaxels, with white corresponding in this case to velocity. Lya: 
flux contours are overplotted. The apertures selected to extract 
the spectra of components Al and A2 (see text) are shown with 
green lines. 




5000 n 5500 



Figure 11. MRC 0140-257: Spectra of Al and A2. The Lya spec- 
tral region is amplified in the right panels. The complete spectrum 
is shown for A2 (left panel) to highlight the detection of CIV and 
Hell (small box). All lines are characterized by a very narrow 
component of FWM<300 km s _1 . The lines are unusually nar- 
row for a high redshift radio galaxy. Notice that the CIV doublet 
is resolved. Flux in units of 10 -16 erg s _1 cm -2 A -1 . 



when both spectra are added. The fit to the line profile in the 
coadded spectrum of Al and A2 is shown in Fig. 12 together 
with the individual components isolated in the fit. The un- 
derlying broad component has FWHM=1200±200 km s -1 . 
The velocity shift between Al and A2 is 120±20 km s"\ 
which is also rather low co mpar ed with typical values in 
HzRG (e.g. McCarthy et al. Il996h . 

CIV and Hell are detected in A2 (Fig. 11, left). Two 
very narrow emission lines are identified at the expected 
wavelengths of the CIV doublet. We have fitted the lines 
with no constraints applied. The resulting doublet consists 
of two emission lines separated by 8.0±1.0 A (consistent 



, A1+A2 














Try r » m ' ,v f 








4380 4400 4420 4440 44t 

X 

Figure 12. MRC 0140-257: Lya spectrum of the coadded spectra 
of Al and A2. The original data are shown with the best fit (thin 
solid line) and the residuals (green). The individual kinematic 
components (blue and red) are shown in the right panels. A strik- 
ingly narrow component of FWHM=270±40 km s _1 is isolated on 
top of an underlying broad components with FWHM =1200 ±200 
km s _1 . Flux in units of 10 -16 erg s _1 cm -2 A -1 . 



within the errors with the theoretical 9.2 A at the redshift 
of the object) and with flux ratio of — 1.0=b0.2, as expected 
in the optically thin case. Taking the errors into account, 
both lines are unresolved (observed FWHM=2.5=b0.6 and 
2.4±0.7 A respectively), with FWHM<100 km s" 1 . This is 
in excellent agreement with the width of the Lya narrow 
component. Hell is also detected and is similarly narrow 
with FWHM=180±60 km s _1 . 

The small FWHM values of the CIV and Hell lines 
confirm that Lya is intrinsically very narrow (rather than 
absorbed) . 

Lya absorption 

The Lya spectral profile does not show clear evidence 
for absorption. This is further supported by the large Lya 
ratios. In A2 Lya/Hell and Lya/CIV are 14.0±1.5 and 
7.4±1 respectively. These values are consistent with stan- 
dard photoionization model predictions for the measured 
CIV/HeII=1.8±0.3 (Villar-Martm et al. 2007). Lya/Hell 
and Lya/CIV are > 6 and >5 respectively for Al and it 
is not possible to say whether Lya is absorbed. The values 
for the integrated spectrum are ^15 and 9±2 respectively 
which for the measured CIV/HeII>1.5 do not imply absorp- 
tion e ither. This is in contradiction with lEales &; Rawlingsl 
(|l996h (see §3.3). 



5 DISCUSSION 

5.1 Main properties of the Lya nebulae 

The three radio galaxies investigated here are associated 
with giant (>60 kpc) Lya nebulae. This was already known 
for MRC 1558-003 and MRC 2025-218, but not for MRC 
0140-257. The total Lya luminosities are -(0.3-3.4) x 10 44 
erg s _1 , within the range of typical values measured for 
HzRG. 

The morphologies are varied. The nebula of one source 
has a centrally peaked, rounded appearance (MRC 1558- 
003). In the other two objects it consists of two spatial 
components. This bimodal morphology could be a conse- 
quence of an obscuring dust structure (e.g. Knopp & Cham- 
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bers 19971 . Reuland et al. l2003h . For MRC 2025-218 this is 
supported by the fact that the conti nuum centroid co i ncides 
with a minimum in the Lya flux. IPentericci" et al. I (|200lh 
also proposed the existence of a dust lane in MRC 0140-257 
to explain the near infrared continuum morphology. 

The alignment between the Lya nebulae and the radio 
structures is remarkable in MRC 2025-217 and MRC 0140- 
257 (<10°). In MRC 1558-003, the nebula is misaligned by 
30±5°. Similar values have been measured for ot her H zRG 
radio galaxies (McCarthy, Spinrad & van Breugel[l995). 

Based on previous, deeper spectroscopic studies and 
this work, we conclude that, although the high surface 
brightness emission tends to be confined within the radio 
lobes, Lya emission is also detected beyond the radio struc- 
tures in the three objects studied here (see also Villar-Martm 
et al. 2003 for other examples). 

In one case (MRC 1558-003), the new data reveal Lya 
emission outside any plausible ionization con e and far from 
the r adio s tructures (see Reuland et al. 120031 Villar-Martm 
et al. 120061 for other examples). Projection effects might be 
responsible, since this is a broad line object and this is our 
favoured explanation (see §5.3). Alternatively, part of the 
Lya emission might be resonantly scattered or powered by 
a mechanism not related t o the active nucleus, such as young 
stars (Villar-Martm et al. l2007h . Unfortunately, the data are 
not deep enough to check whether lines other than Lya are 
emitted in t hese r egions. Cooling radiation (Haiman, Spaans 
& Quataert l2000h is an interesting possibility ofte n disc ussed 
in the subject of radio quiet (e.g. Nilsson et al. 2006) and 
radio loud Lya nebulae. However, the Lya surface brightness 
we measure in regions outside the reach of the ionization 
cones (>10 -17 erg s _1 cm -2 arcsec -2 ) is too high compared 
with the model predictions (see Villar-Martm et al. 120031 for 
a more detailed d iscuss ion on this issue; see also Dijkstra, 
Haiman & Spaans [2Q06h . 

5.2 Interactions between the gas and the radio 
structures 

Interactions between the radio structures (jet and radio 
lobes) and the ambient gas are known to have a profound 
impact on the kinematic properties of the giant nebulae as- 
sociated with many HzRG (§1). 

For the three objects studied here, this is obvious only in 
MRC 2025-218, in the northern region in particular, where 
the emission lines, Lya, Hell and CIV have FWHM>1000 
km s _1 (§4.2). No evidence for interactions is found in the 
southern region, where the radio structures present a sharp 
bend (see Fig. 4). If this has been caused by the collision of 
the radio structures with a dense gaseous regiorjfl, no clear 
signature has remained in the gas kinematic properties. 

Lya is also very broad across the MRC 1558-003 nebula. 
However, FWHMM000 km s" 1 values are measured also in 
distant regions from the radio structures. Since moreover 
Hell is relatively narrow (45 0-650 km s _1 across the nebula, 
see also Villar-Martm et al. l2003T h we cannot discard that 
resonance scattering effects are responsible for broadening 
the Lya spectral profile. 

3 Since this is a broad line object, projection effects might exag- 
gerate the angle of the observed bend of the radio structures 



In MRC 0140-257 the nebular emission is dominated 
by quiescent gas (§4.3), i.e., not perturbed by the radio 
structures. We cannot tell whether the faint Lya under- 
lying broad component (FWHM=1200±200 km s -1 ) is a 
consequence of radio/gas interactions or resonance scatter- 
ing effects. In fact, the most interesting result for this object 
is the strikingly quiescent kinematics revealed by the bulk 
of the emission lines across the nebula (FWHM<300 km 
s _1 and Voffset ^120 km s _1 ). This characteristic is unique 
among HzRG. Similar relaxed kinematics is rather extreme 
also among low z radio galaxies, since only a few show such 
narrow line widths, usually measured in extranuclear regions 
(e.g. Baum, Heckman & van Breugel Q990|) ■ 

5.3 Are the nebulae infalling? 

iHumphrev et al. (2007) proposed that the extended qui- 
escent ionized nebulae associated with numerous powerful 
radio galaxies at different redshifts are in the process of in- 
fall. These authors searched for correlations between several 
radio and optical/UV asymmetries, and found that (i) the 
quiescent ionized gas has its highest redshift on the side 
of the nucleus with the brightest and more polarized radio 
hotspot, and that (ii) on the side where the Lya emission is 
brightest relative to the other emission lines and continuum. 
They concluded that orientation effects, with the quiescent 
gas in infall towards the nucleus, is the most natural scenario 
for explaining these correlations. 

Our study of MRC 2104-242 based on integral field VI- 
MOS data also suggested that the giant quiescent Lya neb- 
ula associated with this object could be in the process of 
infall (Villar-Martm et al. 2006). However, due to the uncer- 
tainty on the spatial gas distribution, a rotational pattern 
could not b e discarded. It was our later work (Humphrey 
et al. 120071 ) on a larger sample what allowed us to discard 
rotation and favour the infall interpretation. 

We investigate next whether the morphological and 
kinematic properties of the nebulae studied here are con- 
sistent with infall. 

MRC 1558-003 

According to IHumphrev et al. I (|2007T ) the giant neb- 
ula associated with MRC 1558-003 is infalling towards the 
center. Here we describe an infall model that can explain 
the observations quite well. An AGN sits in the center of a 
halo of mass M-tot = 5 x 1O 12 M0. A dusty torus surround- 
ing the AGN causes it to photoionize a biconical region of 
the surrounding gas. The opening angle of each cone is as- 
sumed to be cj) — 90°. The cone axis is not perfectly aligned 
with the line-of-sight (L.O.S), but intersects it at an angle of 
~ 20°. Furthermore, the cones are rotated by ~ 10° counter- 
clockwise on the sky. In this scenario, one cone is pointing 
almost directly at the observer, while the other points in 
the opposite direction. Photoionization and subsequent re- 
combination in the cones converts ionizing radiation emitted 
by the AGN into Lya, an d thus results in spati ally extended 
fluorescent Lya emission (|Haiman fe Reesll200lh . In this pic- 
ture, infalling gas in the nearest cone will be emitted with a 
redshift relative to the systemic redshift. 

We found that a model in which the gas density in- 
creases as p oc r -2 , and the gas velocity increases v(r) oc 
r -1 / 2 can reproduce most observed properties of the nebula 
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Figure 13. Observable properties of a model in which an AGN photoionizes a biconical region of the surrounding gas. One cone is 
pointing almost directly at the observer, while the other points in the opposite direction. The gas is collapsing on to the AGN (for model 
details, see Appendix A). If the Lya flux from the cone that points away from the observer is not detected, then this model reproduces 
most observational features (except the Lya FWHM) remarkably well: The left panel shows the surface brightness profile. The contour 
levels are (0.03,0.01.0.2,0.7) times the maximum surface brightness level. The central panel shows the velocity shift of the Lya line with 
respect to the centroid (labels denote the off-set in km s — 1 ). The right panel shows the FWHM in km s — 1 of the Lyo; line as a function 
of position. 



(Fig. 13). In order to match the observations, it is crucial 
that the emission from the furthest cone is not observed 
(see below). For a more detailed description of the model, 
the reader is referred to Appendix A. According to Fig. 13, 
our infall model produces the following features: 

• A rounded and centrally peaked morphology (left 
panel). This is a consequence of the geometry of our model 
and the increase of the recombination rate (and thus of flu- 
orescent Lya emission) towards the AGN. The innermost 
contour level in Fig. 13 encloses a surface brightness that is 
25 times higher than that surrounded by the outermost con- 
tour. These properties are similar to those observed (Fig. 2), 
within a similar nebular size. In order to obtain a more 
symmetric appearance, as shown by the data, the axis cone 
should be closer to the line of sight. This would correspond 
to a quasar. The fact that the AGN continuum and broad 
line emission are detected at rest frame optical wavelengths, 
but are not obvious in the UV rest frame could be explained 
by the existence of dust obscuring the central regions. 

• The largest redshift is observed at the Lya centroid 
(central panel), as observed in the data (Fig. 3, right). Be- 
cause the infall velocity increases towards the AGN, the red- 
shift of the Lya line decreases outwards. Labels in the figure 
denote the mean velocity shift of the emission line in km s _1 
with respect to he centroid, i.e., the position of maximum 
Lya flux. These values are consistent with the measured ve- 
locity offsets. 

• The nebular centroid has the largest redshift w.r.t the 
systemic redshift. It is not possible to determine the sys- 
temic redshift in this object. However, it is interesting to 
note that in MRC 1558-003 Lya is redshifted relative to the 
main rest-frame UV emission lines (CIV, Hell, CHI D both 
in the spatially integrated spectrum (Rottgering et al. 1997) 
and at different spatial positions (Villar-Martm et al. 120031 ). 
In our scenario, the Lya emission from the far cone is neg- 
ligible. The other emission lines might also be fainter from 
the more distant cone because of extinction, but these ef- 
fect would be less pronounced, since they are not affected 
by resonance scattering effects (see also Humphrey et al. 



2007). As a result, the relative contribution from the far 
cone to the flux of these lines is larger than for Lya both 
in the integrated spectrum and at different projected spa- 
tial positions. Since the emission from the more distant cone 
must be blueshifted relative to the near cone because of the 
infall pattern, the lines will be blueshifted relative to Lya, 
which is consistent with the observations. 

To reproduce these two last features, the Lya emission 
detected from the furthest cone must be negligible as the 
expected blueshift of these photons would eliminate the pre- 
dicted increase of the Lya redshift towards the nebular cen- 
troid and the blueshift relative to other emission lines. Al- 
though radiation blueward of the rest-frame Lya frequency 
is subject to absorption in the IGM, the IGM at z — 2.6 
is not opaque enough to completely eliminate the flux from 
the far cone (we used the model of Dijkstra et al. 120071 to 
calculate the impact of the IGM on the Lya line). Alter- 
natively, the flux from the far side could be eliminated by 
a neutral and dusty spatially extended structure that lies 
between the two cones. Several studies suggest the existence 
of su ch structure (e.g. Humphrey et al. 120061 . van Ojik et al. 
1 19971 ). 

Our model has greater difficulty reproducing the ob- 
served Lya FWHM (right panel, Fig. 13). Especially if 
purely radial infall is considered, the model produces a max- 
imum FWHM of -500 km s"\ which is a factor of almost 
3 short of what is observed for Lya. The model FWHM 
shown in Fig. 13 is boosted to 750 km s _1 by convolving 
the spectrum at each pixel with a Gaussian with a standard 
deviation of am — Vdrc/V2 ~ 250 km s _1 (am is the 1-D 
velocity dispersion of the halo). This reflects that the infall 
probably does not occur purely along radial paths. Although 
inconsistent with the Lya measurements, this model repro- 
duces successfully the spatial variation of the Hell FWHM 
observed by (Villar-Martm et al. 2003) along the radio axis 
using long slit spectroscopic data. In this direction, the neb- 
ula shows broader FWHM— 750 km s _1 at the spatial cen- 
troid and —500 km s _1 in the outer parts (Villar-Martm et 
al. 2003), consistent with the models. 
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Resonant scattering of Lya photons has been ignored 
in our model because the gas within the cones is highly ion- 
ized and likely optically thin to Lya (see Appendix A). Al- 
though resonant scattering effects could explain the large 
Lya FWHM values, and its broader profile compared with 
Hell (§4.1), they would contradict oth er observational re- 
sults. [Dijkstra]^aima^^[sj^aa3 ([20061 ) have performed ra- 
diative transfer calculations through neutral collapsing gas 
clouds and found that resonant scattering does not only 
broaden the Lya line, but that energy transfer from the in- 
falling gas to the Lya photons shifts the line to bluer wave- 
lengths. On the contrary, the observations suggest that the 
Lya line is redshifted with respect to other emission lines 
and probably relative to the system's redshift. Neither does 
resonant scattering through optically thick, collapsing, gas 
clouds reproduce the observation that the Lya line is red- 
shifted more towards the center. If resonant scattering does 
not occur, then it remains to be explained why Lya shows 
such large FWHM values and it is much broader than Hell. 

An alternative mechanism that can broaden Lya rela- 
tive to other emission lines and could produce a similar ve- 
locity patt ern as that of infall is outflows (Dijkstra, Haiman 
& Spaans 120061 ). In this scenario, a shell propagates out- 
wards from the center, which is thick enough to be self- 
shielding (with column densities of hydrogen of 10 20 cm -2 ). 
All material is highly ionized, except inside the shell that is 
moving outwards. In this case, Lya photons basically bounce 
off the shell as if it were a mirror, scattering back and forth 
between the two expanding mirrors until they have been 
Doppler boosted far enough for them to be able to simply 
propagate through the shell (each time a Lya photon scatters 
off the shell, it picks up a Doppler shift from the shell). This 
process will broaden Lya, not affecting the other emission 
lines. Detailed modeling would be necessary to test whether 
outflows can reproduce all other observed properties of the 
nebula as infall does. However, since outflows cannot ex- 
plain the radio/op tical asymmetries seen in MRC 1558-003 
(Humphrey et al. l20Q7h . which on the other hand, are suc- 
cessfully explained with infall, here we will not explore the 
outflows scenario further. 

In conclusion, our infall model can explain successfully 
the morphology, size, surface brightness distribution and the 
velocity field of the Lya nebula associated with MRC 1558- 
003. It can also explain why Lya is redshifted relative to 
other emission lines and the FWHM values of the non reso- 
nant Hell line. The infall scenario is also consistent with our 
previous results which imply that the quiescent nebulae in 
many HzRG (and MRC 1558-003 in particular) are infalling 
(Humphrey et al. l200/t ). On the other hand, this model fails 
to explain the large Lya FWHM values. A mechanism which 
might not be connected with the infall process could be re- 
sponsible for the line broadening), while keeping the infall 
pattern and the Hell FWHM intact. 

MRC 2025-218 

It is not possible to investigate the kinematic pattern of 
the quiescent gas, since the data do not allow to isolate its 
emission in the NE component. 

MRC 0140-257 
It is not possible to disentangle how the quiescent gas 



is moving in this object. The relaxed kinematics is not sug- 
gestive of outflows. 

MRC 0140-2 57 do es not follow the trend found by 
Humphrey et al. (2007|) in powerful radio galaxies at dif- 
ferent z, such that the quiescent gas has the highest redshift 
at the side of the brightest and more polarized radio hot 
spot. The authors interpret these results as evidence for in- 
fall. On the contrary, the highest redshift in MRC 0140-257 
is measured in A 1, at the side of the fainter (northern) radio 
hotspot. 

This discrepancy could be due to three different reasons: 
(i) the orientation diagnostics based on the radio properties 
is not valid (ii) the gas is not in infall (iii) our radio-optical 
registration (§2) is incorrect. If the radio core was located on 
Al or A2 (rather than between the two Lya components), 
we would have no information on the redshift on one side of 
the nucleus. 

The interpretation of infall by Humphrey et al. (2007) 
is based on the fact that the radio asymmetries in their 
sample provide a powerful diagnostic tool to constrain the 
orientation of the radio source (and the Lya nebula as a con- 
sequence). One possible explanation for the MRC 0140-257 
discrepancy is that in this particular case the radio proper- 
ties are not reliable indicators of the object orientation. This 
is actually suggested by the fact that the brightest (south- 
ern) hot spot shows the lowest polarization (Carilli et al. 
Il997l ). contrary to what we would expect if both flux and 
polarization asymmetries where mostly determined by ori- 
entation. 

The orientation diagnostics based on the radio proper- 
ties would fail, for instance, when the radio axis is very close 
to the plane of the sky. For MRC 0140-25 7, on t he contrary, 
the detection of a radio core (Carilli et al. Il997h rather sug- 
gests that the angle between the radio axis and the plane of 
the sky is not negligible. 

The diagnostics would also fail if the radio asymmetries 
are a consequence of environmental effects (e.g. McCarthy, 
van Breugel & Kapahi 199l]), rather than orientation. This 
is possibly the case in MRC 0140-257. On one hand, the 
brightest, less polarized hotspot is closest to the nucleus, 
and the Lya emission on that side of the nucleus (adopt- 
ing our present radio-optical registration) is also brightest. 
This is wha t we expect from McCarthy, van Breugel & Ka- 
pahi (1991) environmental scenario. On the other hand, the 
southern (brighter) radio structure consists of two radio hot 
spots (Fig. 8). If this is due to a interactions with the gas 
(e.g. Carilli et al. ll997l ). this process would enhance the ra- 
dio emission and the higher brightness of the southern radio 
structure would therefore be a consequence of environmental 
effects. 

Alternatively, it is possible that the Lya nebula is not 
infalling towards the center. Al and A2 could be two objects 
rotating around a common center of mass (e.g. De Breuck et 
al. 12005). If this is located between Al and A2, we estimate 
a dynamical mass of >3xl0 9 M© for a radius r—lQ kpc 
and a projected rotation velocity of 60 km s -1 (i.e. half 
of the Voffset value between Al and A2). This mass value 
is at least a factor of 100 smaller tha n the stellar masses 
inferr ed for other HzRG (Seymour et al. 120071 Villar-Martm 
et al. 120061 ) or for massive, early- type galaxies at z ~2 in 
deep, IR surveys (e.g. Daddi et al. 2004). MRC 0140-257 
could be an unusually small mass HzRG, a possibility that 
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is also suggested by the small FWHM values measured in 
Al and A2. It would be essential to obtain rest frame H and 
K m agnitudes (e.g. Seymour et al. 120071 . Villar-Martm et al. 
2006) to make a proper estimation of the stellar mass. 

If, on the contrary, this is progenitor of a massive ellipti- 
cal galaxy, as one expects for a powerful HzRG (e.g. McLure 
et al. [1999), the rotation plane of Al and A2 must almost 
coincide with the plane of the sky. 

5.4 Lya absorption 

Lya absorption has definitively been detected in MRC 2025- 
218, spatially extended for >43x32 kpc 2 . The absorber 
could therefore be larger in size than the optical size of the 
host (if already formed) underlying galaxy. CIV is also ab- 
sorbed and, thus, the absorbing gas is highly ionized. As in 
this object, very little kinematic stru cture along the line of 
sight was also found by Wilman et al. (|2005h across the giant 
Lya radio quiet nebula LAB-2 discovered by Steidel et al. 
( 2000|). A more detailed analysis of the absorbers associated 
with MRC 2025-218 will be presented in Humphrey et al. 
(2007, in prep.). 



Lya absorption is definitively detected in MRC 2025- 
218. The absorbing screen has a spatial extension of 43x32 
kpc 2 at least. It is highly ionized (CIV is also absorbed) and 
shows a remarkable velocity coherence (<100 km s _1 ) across 
its spatial extension. 



ACKNOWLEDGMENTS 

Thanks to an anonymous referee for useful comments on 
the paper. We thank Chris Carilli and Laura Pentericci 
for providing the radio maps presented in this paper. The 
work of MV-M and RGD has been supported by the Span- 
ish Ministerio de Educacion y Ciencia and the Junta de 
Andalucfa through the grants AYA2004-02703 and TIC- 
114. SFS thanks the Spanish Plan Nacional de Astronomfa, 
program AYA2005-09413-C02-02 (MEC) and the Plan An- 
daluz de Investigacion of Junta de Andalucfa, research group 
FQM322. AH acknowledges support from a UNAM postdoc- 
toral research fellowship. Thanks to Joel Vernet for produc- 
ing figures 3 and 8. 



6 SUMMARY AND CONCLUSIONS 

The morphological and spectroscopic properties of the giant 
(>60 kpc) Lya nebulae associated with three radio galax- 
ies (MRC 1558-003, MRC 2025-218 and MRC 0140-217) at 
z ~2.5 have been investigated using VIMOS integral field 
spectroscopic data on VLT. The three objects are associ- 
ated with giant (>60 kpc) Lya nebulae. This was already 
known for MRC 1558-003 and MRC 2025-218, but not for 
MRC 0140-257. 

The morphologies are varied. In one source the nebula 
has a centrally peaked, rounded appearance (MRC 1558- 
003), while in two objects it consists of two spatial com- 
ponents (MRC 2025-218, MRC 0140-257). The total Lya 
luminosities are in the range (0.3-3.4) x 10 44 erg s _1 . 

The Lya spectral profile shows strong variation through 
the nebulae, with FWHM values in the range ^400-1500 km 
s _1 and velocity shifts V ff S et ~ 120-600 km s _1 . Kinematic 
disturbance induced by the radio structures plays a clear 
important role in MRC 2025-218. 

Most spectroscopic and morphological properties of the 
giant nebula associated with MRC 1558-003 are successfully 
explained by a scenario such that the giant nebula is col- 
lapsing towards the center in the potential well of a 5xl0 12 
M© halo. This is consistent with our previous conclusion 
that the giant nebulae associated with this and other HzRG 
are m infall (Humphrey et al. 2007). On the other hand, 
this model has problems to reproduce the large Lya FWHM 
values (>1000 km s _1 ), which might be a consequence of a 
mechanism unrelated to the infall process. 

We have discovered a giant (~75 kpc) Lya nebulae as- 
sociated with a radio galaxy at z =2.64. It is very closely 
aligned with the radio axis and it extends well beyond the 
radio structures. It is characterized by a strikingly relaxed 
kinematics (FWHM<300 km s" 1 and V ff S et <120 km s" 1 ), 
unique among HzRG. This object might have an unusually 
small dynamical mass (>3xl0 9 M ) for a HzRG, although 
this needs further investigation with other type of data. 
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gen nuclei at radius r, and we used that the gas is (almost) 
fully ionized. 

In the frame of the cell (which itself is falling toward 
the center), the Lya line is broadened by thermal motions 
of the atoms in the gas, which results in a Voigt profile. The 
bulk motions of the cell then cause a Doppler shift of the 
Lya line. The total spectrum at (a, 5) on the sky, denoted 
by S(a,5, x), is obtained by adding the contributions from 
all cells. To be more precise, the spectrum at (a, S) is given 
by 

S(a,5, x) oc f ds e(a, 8, s) <p( x — \ , (Al) 
J V ^th / 

where e(a,£, s) = nj I (a,5,s) if the cell is inside the cone, 
and otherwise e(a, 5, s) = 0. Furthermore, <j)(x) is the Voigt 
function^. 

The gas density and velocity profiles were assumed to 
be of the forms n#(r) oc r -2 , v(r) oc r _1//2 . The infall veloc- 
ity was assumed to be equal to the circular velocity of the 
halo at the virial radius To prevent both fields from diverg- 
ing, both riH and v were kept constant at radii less than 0.15 
times the virial radius (which corresponds to 0.15 x 160 = 24 
kpc). These fields were chosen to give the best agreement 
with the observations. Note that these parameters were cho- 
sen to produce a good fit to the observations. To enhance 
the FWHM at each location (a, 5) were convolved with a 
Gaussian of width am = ^circ /y/2 ~ 250 km s" 1 . 

Equation IA1I is valid in the limit where the cloud is 
transparent to all Lya photons. In our model the Lya from 
the furthest cone is not observed and we only need to con- 
sider Lya emission from the nearest cone. For these photons 
we find that the optically-thin approximation is very rea- 
sonable: all Lya photons are emitted in the cone where the 
gas is highly ionized. For the cone orientation of our model, 
each Lya photon never leaves the cone on its way to the 
observer. The total optical depth for a given Lya photon 
is reduced even further by velocity gradients along a L.O.S. 
Since all these photons are emitted predominantly redward 
of the Lya resonance, scattering in the IGM is not impor- 
tant. 



APPENDIX A: INFALL MODEL OF MRC 
1558-003 

For each position (a, S) the spectrum was obtained by 
adding the emission from each cell (a, 5, z) inside the cone 
along that line of sight (L.O.S). The emissivity for a cell 
was calculated by assuming photoionization equilibrium. 
This implies the volume emissivity at radius r scales as 
e(r) oc n^(r), where n#(r) is the number density of hydro- 



4 The Voigt function is given by 4>(x) = % (y e _ x y 2 + a 2 , where 

a = A 2 i/4itAi; d = 4.7 x 10 -4 (13 km s~ 1 /v th ). Here A 2 i is the 
Einstein A-coemcient for the transition, vth is the thermal veloc- 
ity of the hydrogen atoms in the gas, given by v t h = ^2h,BT /m p , 
where k,B is the Boltzmann constant, T the gas temperature, and 
m p the proton mass. 



